Spontaneous breaking of rotational symmetry in superconductors* 
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We show that homogeneous superconductors with broken spin/isospin symmetry lower their en- 
ergy via a transition to a novel superconducting state where the Fermi-surfaces are deformed to a 
quasi-ellipsoidal form at zero total momentum of Cooper pairs. In this state, the gain in the conden- 
sation energy of the pairs dominates over the loss in the kinetic energy caused by the lowest order 
(quadrupole) deformation of Fermi-surfaces from the spherically symmetric form. There are two 
energy minima in general, corresponding to the deformations of the Fermi-spheres into either pro- 
late or oblate forms. The phase transition from spherically symmetric state to the superconducting 
state with broken rotational symmetry is of the first order. 



In the original Bardeen-Coopcr-Schricffcr (BCS) the- 
ory of bulk superconductivity, the condensate wave- 
function describes a quantum coherent state which is 
invariant under spatial and time (particle-hole) reversal 
transformations |l| . External perturbations which act on 
discrete quantum variables, like the spin of the fermions, 
break the particle-hole symmetry. A typical example is 
a metallic superconductor in a high magnetic field where 
the Pauli paramagnetism induces an asymmetry in the 
populations of the spin-up and spin-down electrons. The 
superconducting state is quenched via a first order phase 
transition once the splitting in the energy spectrum of 
spin-up and spin-down electrons becomes of the order 
of the pairing gap in the unpolarized state Q- The 
crossover from the BCS to the normal state can be un- 
derstood in terms of the phases space overlap between 
the fermionic states located at the top of their individ- 
ual Fermi-surfaces. The pairing gap is maximal for the 
symmetric (unpolarized) state with perfectly overlapping 
Fermi-surfaces. As these are driven apart by the "po- 
larizing field" the phase-space available for the pairing 
decreases and the gap is successively suppressed. At fi- 
nite temperatures the smearing of the Fermi-surfaces in- 
creases the phase-space overlap and hence the critical 
field at which the superconductivity is quenched. 

The superconducting state sustains larger asymme- 
tries if the translational symmetry is broken. Larkin 
and Ovchinnikov and Fulde and Ferrell (LOFF) argued 
that the crossover from the BCS to the normal state, 
as the Fermi-surfaces are driven apart by the polarizing 
field, occurs via a spatially inhomogeneous superconduct- 
ing phase 0, Q. The Cooper pairs carry a finite total 
momentum in the LOFF phase, i.e., the centers of the 
Fermi-spheres are shifted allowing for a partial phase- 
space overlap. 

This paper suggests an alternative mechanism of 
breaking the symmetry which is based on a deforma- 
tion of the spherical Fermi-surfaces, to the lowest or- 
der, into quasi-ellipsoidal form. If the total momen- 
tum of the Cooper pairs is zero, as we shall assume in 



the following, the deformation spontaneously breaks the 
rotational symmetry. The novel superconducting phase 
maintains its stability due to the dominance of the con- 
densation energy of the Cooper pairs over the loss in 
the kinetic energy of the system caused by the deforma- 
tion of the Fermi-surfaces. We shall assume that there 
is a single axis along which the symmetry is broken, al- 
though more complicated patterns of symmetry breaking 
are possible (simultaneous breaking of the rotational and 
translational symmetries, higher order multipole defor- 
mations of Fermi-surfaces, etc). Note that deformed or 
non-spherical Fermi-surfaces are common for electrons in 
solids; here we treat systems which are homogeneous in 
the normal state, i.e., any deformations of the spheri- 
cal shape of the Fermi-surfaces would correspond to a 
metastable state. 

We start with the BCS-Gorkov equations in energy- 
momentum representation 
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where G 7l g(o;,p) and Fjp(u),p) refer to the full normal 
and anomalous retarded propagators, A a7 (a;,p) is the 
anomalous self-energy, and the diagonal matrix elements 
of the first matrix correspond to the inverse of free-single 
particle propagators; the Greek indeces a, (3... = 1,2 
label the two different species, and w and p refer to the 
particle energy and momentum. (Note that the center-of- 
mass momentum of particles is zero). Suppose that the 
rotational symmetry is broken by a deformation of the 
Fermi-surfaces from spherical form. The quasiparticle 
spectrum of species a can be parameterized, then, as 
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where /i Q are the chemical potentials of the particles in 
the undeformed state, 9 is the angle between the parti- 
cle momentum p and the axis of symmetry breaking; the 
deformation of the Fermi-sphere in Eq. (J2J is truncated 
at the lowest order non-trivial axisymmetric deformation, 
which is described by the I — 2, m — term of the Legen- 
dre polynomials associated with an expansion in spherical 
harmonics. The constant energy surfaces of quasiparticle 
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excitations defined by Eq. © represent quasi-ellipsoids 
of revolution with an ellipticity e a . For e a — the spec- 
trum (j2J is the true eigenstate of the unpaired, homoge- 
neous system in the absence of external fields. (Note that 
the deformation described by Eq. does not need to 
conserve the volume of the Fermi-sphere, as we impose 
a self-consistency condition for the total density of the 
system, see Eq. (0 below). 

In the following we shall neglect the possible pairing 
among the same species (A aa — 0) so that only the off- 
diagonal elements of the anomalous self-energy matrix 
are non-zero. The quasiparticle excitation spectrum in 
the superconducting phase is determined in the standard 
fashion by finding the poles of the propagators in Eq. 
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where the symmetric and anti-symmetric parts of the 
spectrum (which are even and odd with respect to the 
time-reversal symmetry) are defined as E$,a — (Ei ± 
E 2 )/2. The solution of Eq. Q can be written in terms 
of the eigenstates © as 
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where the Bogolyubov amplitudes are 

u 2 = -+ Es v 2 = - ^~ 
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The mean-field approximation to the anomalous self- 
energy yields the gap equation 
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where V(p, p') is the bare interaction 01 > /( w ) — 
[exp(/3o;) + 1] 1 is the Fermi distribution function and /3 
is the inverse temperature. The lo integration is straight- 
forward in the quasiparticle approximation, since the fre- 
quency dependence of the propagator is constrained by 
the on-shell condition. For S-w&ve interactions the po- 
tential depends only on the absolute magnitude of the 
quasiparticle momenta. In this case the gap equation 
simplifies to 
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Note that the deformation of the Fermi-surfaces enters 
the gap equation as a parameter which is determined by 
the minimum of the ground state energy of the supercon- 
ducting phase. For strongly coupled superconductors the 



gap equation (J8J is supplemented by the normalization 
condition for the net density p = p\ + p 2 at a constant 
temperature. The densities of the species are given by 
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where the summation is over the discrete quantum vari- 
ables. The second equality follows in the quasiparticle 
approximation. 

Next we turn to the thermodynamic properties of the 
superconducting state. At a fixed density and temper- 
ature the relevant thermodynamic potential is the free 
energy: 



F\ PtP = U-0- 1 S, 
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where U is the internal energy and S is the entropy. In 
the mean-field approximation the entropy is given by the 
expression 

S = -ks^J ^{/(w 1 )ln/( Wl ) + /> 1 )ln/> 1 ) 

+ /(w 2 )ln/(w 2 ) + /(w 2 )ln/(w 2 )}, (11) 

where f{u)±) = [1 — f(u>±)], ks is the Boltzmann con- 
stant. The internal energy, defined as the grand canoni- 
cal statistical average of the Hamiltonian, is 
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The first term in Eq. I|12|l is the kinetic energy while the 
second term includes the mean field interaction among 
the particles in the condensate. The true ground state 
of the system minimizes the free energy difference 8F\ p .p 
between the superconducting and normal states (the free 
energy in the normal state follows from Eqs. (|llfl and 
when A = 0). 
The deformations of the Fermi-spheres can be de- 
scribed in terms of the "conformal deformation" e = (ei + 
e 2 )/2 and the "relative deformation" Se = (ei — e 2 )/2. 
Then, the symmetric and anti-symmetric parts of the 
energy spectrum can be written as 
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where p = (p\ + /i2)/2, Sfi — (fix — Hzj/I (here we ignore 
the difference in the masses of the spin/isospin up and 
down quasi-particles) . Equations (jBJl, © an d ifTH)) form 
a closed system, which determines the pairing gap and 
the ground state energy of a superconductor for constant 
density asymmetry a = (pi — p2)/{pi + Pt). The values 
of the deformation parameters Se and e are obtained by 
requiring that the free-energy attains its minimum. Note 
that in the weak coupling limit Eqs. © and © decouple 
and one may solve for A as a function of Sp instead of 
a. Apart from the fact that a, rather than Sp, is the 
measurable quantity, there is an additional reason for 
solving the full set of equations. The gap equation alone 
is symmetric under exchange Ea — * —Ea, which implies 
that the solutions are symmetric under the simultaneous 
change of the signs of Sp and Se. Eq. ©, however, does 
not have this symmetry and the solutions are distinct 
under the sign transformation above. 

As a specific example, which illustrates the solutions 
above, we consider isospin-singlet (neutron-protonVpair- 
ing in nuclear matter in the 3 Si- 3 Di channel 0, 0, 0- 
The gap in the isospin symmetric case is Aoo = 12 MeV if 
we use as the bare interaction the Argonne potential and 
ignore the renormalization of the mass of the particles in 
the normal state due to interactions. The modification 
of the particle self-energy in nuclear medium (see for a 
review |8() affect the absolute magnitude of the gap and 
rescale its dependence on the parameters. Clearly, with 
these approximations, our model is schematic, however 
we do not expect qualitative changes when renormaliza- 
tion of the interaction and the bare mass are included. 

The BCS solutions for the n-p pairing has been stud- 
ied for the homogeneous (translationally and rotation- 
ally) invariant state under isospin asymmetric condi- 
tions jij 0, 0, ^| and the inhomogeneous state with 
broken translational symmetry (l3llT4| (the nuclear ana- 
log of the LOFF phase. (The flavor asymmetric < qq >- 
condensate in high density QCD is another example of 
strongly coupled superconductor where the breaking of 
translational symmetry plays a role |15|). Consistent 
with the assumption A aa = above we ignore the n-p 
pairing in the Sq channel fl6| . 

Fig. 1 shows the pairing gap as a function of the den- 
sity asymmetry and the relative deformation Se for van- 
ishing conformal deformation (e = 0). Here and below we 
assume a density p = 0.16 fm -3 and a temperature T = 3 
MeV. The gap parameter is normalized to its value Aoo 
in the isospin symmetric rotationally/translationally in- 
variant state. Although a changes in the interval [— 1; 1] 
in general, the symmetry of the equations with respect 
to the indexes labeling the species reduces the range of a 
to [0; 1]. Assuming neutron excess implies that the Greek 
indeces equal 1 refer to neutrons and while those equal 
2 refer to protons. The relative deformation obviously is 
not bounded and can assume both positive and negative 
values. For positive values of Se, which imply an oblate 
deformation for the Fermi-surface of neutrons and a pro- 
late deformation for the Fermi-surface of protons, the so- 




FIG. 1: The pairing gap as a function of the density asymme- 
try a and the relative deformation Se. The gap is normalized 
to its value for Sot — = Se. 
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FIG. 2: A projection of the Fermi-surfaces on a plane parallel 
to the axis of the symmetry breaking. The concentric cir- 
cles correspond to the two populations of spin/isospin-up and 
down fermions in spherically symmetric state (Se = 0), while 
the deformed figures correspond to the state with relative de- 
formation Se = 0.64. The density asymmetry is a — 0.35. 



lutions for the gap equation show the following features. 
For arbitrary constant 5e the gap is maximal at a = 
and is suppressed as the asymmetry is increased. For 
constant a, dA/dSe = corresponds to a maximum at 
Se ^= in the large a limit. The position of the maximum 
is independent of a and is located around Se = 0.5 in our 
model; this value also corresponds to the critical asym- 
metry a c at which the superconducting state vanishes. 
Note that for a around a c the gap exists only in the de- 



4 



formed state. For a = 0, Eqs. (JSJ-JTOJ are symmetric 
under exchange of the sign of deformation and so is the 
gap function. In particular, for a = 0, the critical defor- 
mation for positive and negative deformations coincide. 
For finite a the dependence of the gap on the relative 
deformation depends on the sign of Se. In contrast to the 
positive range of Se, where the maximum value of the gap 
is attained at constant Se, for negative Se the maximum 
increases as a function of the deformation and saturates 
around Se ~ 1. Quite generally, to maintain the maximal 
phase space overlap, the system prefers to keep the sign 
of 5/j, opposite to that of Se. Fig. 2 shows a 2-dimensional 
projection of a configuration of deformed Fermi-surfaces 
for Se — 0.64 which minimizes the free-energy for fixed 
a = 0.35. 




FIG. 3: The difference in the free-energies of the supercon- 
ducting and normal states SF as a function of the density 
asymmetry a and the relative deformation Se. The free-energy 
is normalized to its value for 8a — Se = 0. 

The difference in the free-energies of the superconduct- 
ing and normal states SF is shown in Fig. 3. Owing to 



the symmetries of the underlying equations, SF is sym- 
metric with respect to the sign change of Se when a = 0. 
For finite a's, a minor departure from the rotational in- 
variant state leads to a decrease in SF which develops 
two minima for either sign of Se. This behavior can be 
traced back to the increase of the potential energy with 
increasing gap (cf. Fig. 1). Note that, although there are 
non-trivial solutions to the gap equation in the large a 
limit for Se —* — 1, these solutions do not lower the energy 
of the system. The superconducting phase becomes un- 
stable for a > 0.4 due to the increase in the kinetic energy 
caused by the large deformation of the Fermi-surfaces, so 
that SF is a nearly even function of Se. An inspection 
of the latent heat associated with the phase transition 
at finite temperatures shows that this quantity does not 
vanish at the crossover from the spherically symmetric 
to the deformed superconducting state. Consequently, 
the phase transition associated with the breaking of the 
rotational symmetry is of the first order. 

To summarize, this paper suggests a novel mechanism 
of symmetry breaking in superconducting system with 
particle-hole asymmetry. The lowest order (quadrupole) 
deformation of the Fermi-surfaces (at zero total momen- 
tum of the Cooper pairs) increases the phase-space over- 
lap between the Fermi-surfaces of paired quasiparticlcs, 
which is otherwise depleted by the asymmetry in the par- 
ticle/hole populations. As a result, the free-energy de- 
velops minima for finite deformations, since the gain in 
the (negative) pairing potential energy dominates the in- 
crease in the kinetic energy caused by the deformation. 
Since the deformed ground state spontaneously breaks 
the rotational symmetry, the dynamic properties of the 
superconducting state with deformed Fermi-surfaces such 
as the sound attenuation, the infrared absorption or the 
Meissner effect will be anisotropic. The results above do 
not depend on the nature of fields inducing the asymme- 
try in the fermion populations nor on the nature of the 
pairing forces and should be applicable to a wide range 
of fermionic systems. 
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